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Abstract: An interesting feature of the next-to-minimal supersymmetric standard model 
(NMSSM) is that one or more Higgs bosons may be comparably light (M^ < Mz) without 
being in conflict with current experimental bounds. Due to a large singlet component, their 
direct production in standard channels at the Large Hadron Collider (LHC) is suppressed. 
We demonstrate that there are good prospects for observing such a light Higgs boson 
in decays of heavy neutralinos and charginos. We consider an example scenario with 
20 GeV < Mjj-i < Mz and show that a large fraction of the cascade decays of gluinos and 
squarks involves the production of at least one Higgs boson. Performing a Monte Carlo 
analysis at the level of fast detector simulation, it is demonstrated how the Higgs signal 
can be separated from the main backgrounds, giving access to the Yukawa coupling of 
the Higgs to bottom quarks. Analyzing the resulting bb mass spectrum could provide an 
opportunity for light Higgs boson discovery already with 5 fb" 1 of LHC data at 7 TeV. 
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1 Introduction 

The precise nature of the Higgs mechanism thought to be responsible for electroweak 
symmetry breaking remains unknown. To discover and study the properties of one or 
more Higgs bosons is therefore a challenge — and one of the major objectives — for the 
experiments at the running Large Hadron Collider (LHC). Most of the Higgs boson search 
strategies to date are designed to probe the Higgs sector of the Standard Model (SM), or 
its minimal supersymmetric extension (MSSM) [1, 2]. 

There are several theoretically appealing arguments for weak-scale supersymmetry to 
be realized in nature: it solves the hierarchy problem of the SM Higgs mass, it enables gauge 
coupling unification, and with R-parity conservation it also provides a natural dark matter 
candidate. On the other hand, the realization of weak-scale supersymmetry in terms of the 
MSSM is not free of theoretical problems, such as the scale for the bilinear //-parameter 
entering the MSSM superpotential with positive mass dimension. This parameter has 
no natural values besides Mqtjt or zero, while at the same time it must be close to the 
electroweak scale for a phenomenologically acceptable theory. To solve this problem in an 
elegant way the MSSM can be extended by a complex scalar singlet, giving the so-called 
next-to- minimal supersymmetric model (NMSSM). In this model an effective /i-term of 
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the right size can be generated dynamically from supersymmetry-breaking operators. For 
a general introduction to the NMSSM we refer to the recent reviews [3, 4]. 

The NMSSM is characterized by an enlarged Higgs and neutralino sector as compared 
to the MSSM, giving rise in particular to a richer Higgs phenomenology While it is well 
known that in certain scenarios of the MSSM with complex parameters a light Higgs, with 
mass much below that of the Z boson, is unexcluded by the searches at LEP [5] and the 
Tevatron [6] (see [7, 8] for recent reevaluations with improved theoretical predictions), such 
a scenario can occur even more generically in the NMSSM. In order to be compatible with 
the limits from the LEP Higgs searches, in particular the couplings of such a light Higgs 
state to gauge bosons must be heavily suppressed. As a consequence of the presence of a 
Higgs singlet, in the NMSSM such a situation happens whenever a light Higgs state has a 
sufficiently large singlet component. 

The search for a heavier Higgs state with SM-like (or only moderately suppressed) 
couplings to gauge bosons is complicated in such a scenario by the fact that often the 
decay of this heavier Higgs state into a pair of lighter Higgses is kinematically open, giving 
rise to unusual decay signatures and to a large suppression of the standard search channels 
for a SM-like Higgs. It should be noted in this context that the observation of a decay 
of a heavier Higgs into a pair of lighter Higgses would provide an opportunity for gaining 
experimental access to triple-Higgs couplings, which are a crucial ingredient of electroweak 
symmetry breaking via the Higgs mechanism. 

While in the NMSSM the case of a very light pseudo-scalar, Ma 1 < 2m&, has found 
considerable attention in the literature, in particular in the context of "ideal" Higgs sce- 
narios [9-14], we will focus in the following on scenarios with a light CP-even boson with 
20 GeV < M Hl < M z . Within the MSSM the best known example of a light Higgs 
that is unexcluded by the present search limits is the "hole" in the coverage of the CPX 
benchmark scenario [15] for Mh 1 ~ 45 GeV and moderate values of tan/3 [5], see [8] for 
a detailed discussion of the dependence of the unexcluded parameter region on the choice 
of the various MSSM parameters. It will be difficult to cover this parameter region with 
the standard search channels at the LHC [16-18]. Various other (non-standard) search 
channels have been proposed which may provide additional sensitivity in the quest to close 
this "CPX hole" [16, 19-26]. 

In our analysis within the NMSSM we will investigate the prospects for the production 
of light Higgs bosons in cascade decays of heavy SUSY particles at the LHC. Such an 
analysis, where Higgs bosons are produced in association with — or in decays of — other 
states of new physics, is necessarily more model-dependent than the Higgs search in SM- 
like channels. On the other hand, investigating Higgs physics in conjunction with the 
production of other states of new physics offers additional experimental opportunities and 
may also be more realistic, in the sense that in order to extract a Higgs signal backgrounds 
both from SM-type and new physics processes have to be considered. In the case of the 
MSSM with real parameters, for a Higgs with a mass above the LEP limit for a SM Higgs 
of 114.4 GeV [27], a detailed experimental study for Higgs boson production in a SUSY 
cascade has been carried out by the CMS Collaboration [28], involving a full detector 
simulation and event reconstruction. These results, obtained for the benchmark point 
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LM5, cannot be directly translated to the case of searches for a Higgs boson with mass 
far below Mz, since in the latter case the b jets resulting from the Higgs decay tend to be 
softer. Further phenomenological analyses of Higgs production in SUSY cascades in the 
MSSM with real parameters (and Higgs masses above the LEP limits) have been carried 
out in [29-31], with recent developments focusing on jet substructure techniques to identify 
highly boosted Higgs bosons and enhance the discovery significance [32, 33]. The case of 
a lower mass Higgs has been considered in [23], and it has been pointed out that in the 
CPX scenario there is a significant rate for producing a light MSSM Higgs boson in SUSY 
cascades, but no simulation of signal and background events was performed. The potential 
importance of SUSY cascades to establish a signal for a light CP-odd Higgs in the NMSSM 
has been pointed out in [34]. 

We generalize and extend the investigations carried out in [23, 34] by calculating the 
sparticle decay modes in a general NMSSM setting and performing a Monte Carlo simula- 
tion of the signal and the dominant background to the level of fast detector simulation. A 
simple cut-based analysis is performed to demonstrate that signal and background can be 
resolved in the 66+jets channel. The observation of the Higgs decay in the bb final state 
would be of interest also as a direct manifestation of the Higgs Yukawa coupling. 

The outline of our paper is as follows: the next section begins with a brief recapitulation 
of the NMSSM, presenting the scenario with a light CP-even Higgs boson in some detail. 
In section 3, we describe the production of squarks and gluinos at the LHC and their 
eventual decay into Higgs bosons through electroweak cascades involving neutralinos and 
charginos. Section 4 describes a phenomenological Monte Carlo analysis of these cascade 
processes and contains the main results of this work in terms of kinematic distributions 
demonstrating the separation of signal and background. The conclusions are presented in 
section 5. 

2 The Next-to-Minimal Supersymmetric Standard Model 

In this section we review briefly the elements of the NMSSM which differ from the MSSM. 
Our conventions for the other sectors — that remain unchanged when going to the NMSSM 
— follow those of [35, 36]. 

2.1 The Higgs sector 

The ^-symmetric version of the NMSSM is given by the scale-invariant superpotential 
^nmssm = Yu q l . r Jjc l + y d Q L ■ H d D c L + Y e L L ■ H d E c L + \SH U ■ H d + ^kS 3 , (2.1) 

where <3? denotes a chiral superfield with scalar component The complex gauge singlet S 
is a new addition with respect to the MSSM. To have a complete phenomenological model 
the soft SUSY-breaking terms must also be specified. These are extended by couplings of 
the singlet field, giving new contributions to the scalar potential 

y NMSSM = y MSSM + ^2 + . + ^ kAk S 3 . (2.2) 
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The NMSSM Higgs potential, which is derived from the usual .F-terms, D-terms and the 
soft-breaking potential given by Equation (2.2), allows for a minimum where the singlet 
develops a vacuum expectation value (vev) v s = (S). This induces an effective bilinear 
term A (S) H u ■ H^, thus providing a dynamical explanation for the fi parameter of the 
MSSM in terms of fj, e g = \v s . 

Electroweak symmetry breaking (EWSB) proceeds similarly to the MSSM, and the 
two Higgs doublets are expanded around the potential minimum according to 

/^+^-^>v Hu= ( « y (2. 3) 



Equivalently, the singlet field has an expansion 



S = v s + -= (4> s + ia s ) . (2.4) 

Using the minimization conditions for the potential, the scalar mass parameters w,# u , m H d ' 
and m| can be traded for 

M 2 =g 2 v 2 = \(g 2 +g 2 )v 2 , 
tan/3 = v u /v d , 

fi ef[ = \V S , 

where the doublet vevs fulfill v 2 = v\ + v 2 d = (174 GeV) 2 . Assuming a CP-invariant Higgs 
sector, all parameters are taken to be real. The number of parameters is increased from 
the MSSM. In addition to Ma (or M H ±), and tan/3, the values of A, k, and A K can be 
chosen as free parameters. 

After EWSB, the addition of a complex scalar field gives rise to additional particles 
in the NMSSM spectrum with respect to the MSSM: two additional Higgs bosons (one 
of which is CP-even, the other CP-odd) and their fermionic partner, the singlino. The 
elements of the tree- level mass matrix M 2 H for the CP-even Higgs bosons are given in the 
basis ((f>d, (f>u, (f>s) by 

(M 2 H ) n = Mf cos 2 P + Bfi cS tan /3, 
{M 2 H ) 22 = Mf sin 2 $ + B(j, eS cot 0, 
A ^ A- 

(M 2 H ) 33 = cos/3sin/3 + kv s (A k + 4kv s ) , 

(M 2 H ) 12 = (2AV - Mf) cos sin (3 - B^l cS , 
{M- 2 h) 13 = [2/Ucfj cos P — (B + kv s ) sin /3] , 
(•Mh) 2 3 = ^ v I 2 ^efT sin/3 - (B + nv s ) cos P] , 

where B = A\ + kv s . This matrix is diagonalized by a real 3x3 matrix with elements Sij , 
such that the Higgs mass eigenstates Hi are given by Hi = Sijcfij. For the CP-odd states, 
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the mass matrix elements in the basis (a^, a u ,a s ) can be written 
(M 2 A ) U = fi/iefftan/3, 



{M 2 A ) 22 
(M 2 A ) 12 

13 



(B + 3kv s ) cos /3 sin /3 — 3K^4 K t; s , 



(2.6) 



B/j, cS cot /3, 
A^ 

B/J-eff, 

Xv (B — 3kv s ) sin /3, 
(.B — 3kv s ) cos /3. 

Similarly to the CP-even case, the massive eigenstates Ai can be written using a mixing 
matrix as Ai = Pij&j- One degree of freedom is massless and corresponds to the neutral 
Goldstone boson providing the longitudinal degree of freedom to the Z boson. For some 
purposes it can be convenient to introduce the CP-odd mass parameter 
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Mi 



C/3S/3 



(2.7) 



which corresponds to the mass of the CP-odd Higgs boson in the MSSM limit of the 
NMSSM. 1 

No additional charged scalar is introduced in the NMSSM, but the relation of the 
physical charged Higgs boson mass to the CP-odd mass parameter gets modified. At 
tree-level the charged Higgs mass is now given by 

r2 



M 



aV. 



(2.8) 



2.2 The neutralino sector 



As a result of introducing the new singlet superfield S, the NMSSM comes with an addi- 
tional fermion partner of the complex scalar S, the singlino. The singlino mixes with the 
existing four neutralinos of the MSSM. The resulting 5x5 mass matrix derives from the 
bilinear terms 

C 



(2.9) 



and in the basis (— LB, — iW, Hj, H u , S) it is given by 
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(2.10) 



The upper left 4x4 submatrix is identical to the neutralino mass matrix in the MSSM. 
The neutralino masses can be diagonalized by a single unitary matrix N such that 

D = diag(m^o) = N*M^ a N^ (2.11) 

0, k — y 0, while keeping the ratio k/A and all dimensionful 



lr rhe MSSM limit is obtained by taking A 
parameters fixed. 
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is real and positive with the neutralino mass eigenvalues in ascending order. Alternatively 
one can use a real mixing matrix N, and allow D to have negative elements. In this case 
the physical neutralino masses are given by Im^ol and the neutralino couplings incorporate 
the additional phase shift on the neutralino fields. 

2.3 The squark sector 

We adopt a universal value MgusY for the soft SUSY-breaking scalar mass parameters. 
This means that, for each squark pair qi, q~R of a given flavour, the mass matrix attains 
the form 

M 2 = ( M£ VSY + m\ + M 2 Z cos 2/?(lf - Q q s 2 w ) m q X q \ 

q \ rn q X q M| USY +m 2 g + M 2 Z cos 2PQ q s% r J ' 

(2.12) 

Here m q is the mass of the corresponding quark, 1^ the third component of the weak isospin, 
and Q q the electric charge quantum number. For the weak mixing angle we introduce the 
short-hand notations sw = s'm8w and cy/ = cos 0w- The off-diagonal elements of M"~ are 
related to the soft trilinear couplings A q as X q = A q — fi c g cot [3 for up-type squarks, and 
X q = A q — ^ c fjtan/3 for the case of down-type squarks, respectively. The mass eigenstates 
Q.2) are obtained by a diagonalization of the mass matrix. A generic squark mass will 
be denoted M q below. 

2.4 Scenarios with light Higgs bosons 

As mentioned above, we will focus in the following on the case where the lightest CP-even 
Higgs boson of the NMSSM, H\, has a mass much below Mz- The fact that such a light 
Higgs, possessing a heavily suppressed coupling to gauge bosons as compared to the Higgs 
boson of the SM, may be unexcluded by the current search limits is known already from 
the case of the MSSM with complex parameters [5, 7, 8]. In the NMSSM such a situation 
happens more generically, in particular also for the case where the SUSY parameters are 
real. If the mass eigenstate H\ has a large component of the singlet interaction state cj> s , 
its couplings to gauge bosons (and also to quarks) will be correspondingly suppressed. We 
will investigate the prospects for detecting such a light Higgs state through its production 
in SUSY cascades. 

In the numerical analysis, we shall use a scenario derived from the "P4" benchmark 
point defined in [37]. This benchmark can be realized in models with non-universal Higgs 
mass parameters 7^ mH d ) at the scale of grand unification, and it is compatible with 

the data on the cold dark matter density. As originally defined, the P4 benchmark contains 
a very light CP-even Higgs boson (Mjjx = 32.3 GeV). In order to explore the full range 
Mjji < Mz, we slightly modify the scenario to allow changing the value of Mn x , with the 
remaining phenomenology essentially unchanged. To this end we set A = 0.6 and allow 
A K to take on values in the range GeV < A K < 300 GeV. 2 The soft SUSY-breaking 
parameters are defined directly at the SUSY-breaking scale, allowing us to consider a more 

2 The original P4 benchmark point is recovered for A = 0.53 and A K = 220 GeV. 
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Table 1. Values for the NMSSM input parameters at the SUSY-brcaking scale in the modified P4 
scenario. 
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Figure 1. Higgs mass spectrum for the modified P4 scenario with Msusy = 750 GeV as a function 
of the free parameter A H . For .4,,- > 250 GeV (dashed) the clcctrowcak symmetry remains unbroken 
in the global minimum. 

general spectrum for the remaining (non-Higgs) sectors of the theory. Values for the tree- 
level parameters in the Higgs sector and the soft SUSY-breaking parameters in the modified 
P4 scenario are specified in table 1. The two values given for MgusY will be used later in 
the phenomenological analysis, while the values quoted in this and the next section (unless 
otherwise stated) have been evaluated for Msusy = 750 GeV. 

The NMSSM Higgs masses are subject to sizable corrections beyond leading order 
[38-43]. In order to incorporate the most accurate predictions currently available [44], 
NMSSMTools 2.3.5 [45-47] is used to compute the Higgs spectrum. The resulting Higgs 
masses in the modified P4 scenario are shown in figure 1 as a function of the free parameter 
A K . In the region with A K > 250 GeV the global minimum of the Higgs potential does 
not break the electroweak symmetry; hence these values will not be considered. The 
masses of two Higgs bosons show a dependence on A K : the lightest CP-even Higgs H\ 
(Mh 1 varying from about 110 GeV to 20 GeV), and the lightest CP-odd Higgs A\ (with 
Ma x going from about 90 GeV to 200 GeV). H2 is always SM-like and has a mass of 
about M H2 = 115 GeV (118 GeV for Msusy = 1 TeV). For all Higgs masses in this 
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Figure 2. Branching ratios of iTi (upper left), H2 (upper right), and A\ (lower) in the modified 
P4 scenario with M S usy = 750 GeV. 



plot the NMSSM scenario is compatible with the direct limits from Higgs searches. The 
light CP-even Higgs (Mhj <C Mz) is allowed due to a large singlet component, with |Si3| 2 
ranging from 0.9 for A K = GeV to IS13I 2 > 0.99 for A K = 250 GeV. As a consequence, 
the couplings of H\ to vector bosons are heavily suppressed, so that the cross section for 
production through Higgsstrahlung drops below the LEP limit. The pair production of 
A\H\ is even further suppressed by the large singlet fractions of both Hi and A±, while 
production of H2A1 and H2Z are beyond the kinematic reach of LEP. The full mass ranges 
shown in the figure are also compatible with the constraints from i?-physics implemented 
in NMSSMTools 2.3.5 [45-47], as expected when the charged Higgs boson is heavy [48, 49]. 
The precise values obtained here for the heavy Higgs masses are M#± ~ 563 GeV, and 
Mh 3 — Ma 2 — 572 GeV. None of the heavy Higgs bosons will play any role in the following. 
With the negative sign for the effective /i parameter, this model cannot be used to explain 
the observed deviation in the anomalous magnetic moment of the muon (see e.g. [50] for 
a review) . However — since the considered value of tan /3 is rather low — the predicted 
value for (g — 2) M at least stays close to that in the SM. 

The branching ratios of the three lightest Higgs states, Hi, H2, and A\, are given 
in figure 2. As can be seen from this figure, the light singlet Hi decays preferentially into 
bb, with BR(i7i — > bb) ~ 90% over the full mass range. The subdominant decay into 
tt basically saturates the Hi width. For lower values of A K — where Mh 1 ^ 90 GeV 
- we note a similar enhancement of BR(i/i — > 77) compared to a SM Higgs with the 
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same mass as recently discussed in [51]. The H2 has a more complicated decay pattern, in 
particular for low A K where H2 —> bb dominates and several competing modes (H? — > tt, 
gg, WW) each have a branching fraction around 10%. In this region H2 is SM-like, and the 
same search strategies as devised for the SM Higgs (and the lightest MSSM Higgs boson 
in the decoupling limit) should apply. This situation changes radically when the channel 
H2 — > H\H\ opens. When this is the case, the H2 — > H\H\ mode becomes completely 
dominant. Finally, the lightest CP-odd Higgs A± decays predominantly into bb, with a 
large fraction going into the mode A± — > H\Z when kinematically accessible. 

In the neutralino sector the mass spectrum is independent of A K (and Mstjsy) at lowest 
order, cf. equation (2.10), and therefore remains fixed at: 

M*o = 97.6 GeV, M#> = 227 GeV, M^o = 228 GeV 

Xi Xi X-j, 

M^o = 304 GeV, M^o = 616 GeV. 

There is a clear hierarchy in the mass parameters, which leads to a small mixing between 
the neutralinos. The heaviest neutralino is almost exclusively wino, and Xa is mostly 
bino. The intermediate mass states yQ-> an d X3 are predominantly Higgsino, while the 
lightest neutralino Xi is the singlino. The lightest neutralino is also the overall lightest 
supersymmetric particle (LSP) in these scenarios and thereby a candidate for cold dark 
matter. 

3 Higgs production in the light Hi scenario 
3.1 Standard channels 

The rate for direct production of a light singlet H\ in gluon fusion, gg — > Hi, is proportional 
to its reduced (squared) coupling to quarks. Compared to a SM Higgs boson with the same 
mass, the dominant top loop contribution contains the additional factor |5i2| 2 / sin 2 /3. The 
size of IS12I 2 is limited from above by IS12I 2 < 1 — | -S'is | 2 , where £13 is the singlet component. 
For Mh 1 <C Mz, where | -S13 1 — > 1, the rate for this process gets heavily suppressed. The 
cross section for Hi in weak boson fusion, involving the coupling of Hi to gauge bosons, 
is similarly suppressed. In scenarios where Mh 1 > Mz (corresponding to the mass range 
below the LEP limit on a SM-like Higgs which is unexcluded in the MSSM with real 
parameters) the suppression of gg — > Hi can be overcome by an increased branching ratio 
for Hi -> 77 [51]. 

For A K > 200 GeV in the modified P4 scenario, Hi is light enough to be produced 
through the decay of the SM-like H2 — > H\H\, which can be dominant, see figure 2. The 
production of H2 in standard channels is not suppressed. The resulting two-step decay 
chain leads to "unusual" final states for H2: 46 (about 82% of all decays), 262r (17%), and 
4r (0.6%). These final states make it difficult to establish a Higgs signal, as it has been 
demonstrated, for instance, by the numerous attempts [52-56] to establish a "no-loose" 
theorem for NMSSM Higgs searches when decays of the SM-like Higgs into lighter Higgses 
are open. 
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Another possibility to produce H\ in Higgs decays would be through the decay A\ — > 
H\Z. However, the singlet nature of A\ in the modified P4 scenario leads to a suppression of 
A\ production similar to that for Hi, and this mode is therefore not likely to be accessible. 

The direct production of the heavy Higgs bosons H 3 , A2, and H^ is in principle not 
suppressed with respect to the MSSM case, but at a mass close to 600 GeV and low tan j3 
the observation of those states at the LHC will be difficult even at high luminosity. A large 
fraction of the heavy Higgs bosons in this scenario will decay into lighter Higgs bosons, 
neutralinos and charginos. A detailed investigation of these channels could possibly be of 
interest for a study assuming a very high luminosity at 14 TeV, but is beyond the scope of 
the present paper. 

In summary, it will be problematic to produce and reconstruct the light H\ in any of 
the standard channels proposed for Higgs production at the LHC. We shall focus instead 
on the possibility to produce H\ in the decays of super symmetric particles. 

3.2 SUSY cascades 

As discussed in the previous section, inclusive production of the heavier state H 2 with 
subsequent decay H2 —> H\H\ may be difficult to observe at the LHC. However, the 
related process where a heavier neutralino decays into a lighter neutralino and a Higgs 
boson (and the corresponding mode of the decay of the heavier chargino) may offer better 
prospects. In fact, a light Higgs boson in the mass range below Mz may occur in a 
large fraction of cascade decays of heavier SUSY particles that are produced via strong 
interaction processes. The hard scale associated with the sparticle production can lead to 
event signatures which are more clearly separable from the SM backgrounds than those of 
inclusive Higgs production followed by a decay into a pair of H\ states. The processes of 
interest are 

Xi^XjHk, Xi^XjAk, (3.1) 
xf^xtH k , xf^xtA k , (3.2) 

where H k (A k ) denotes any of the CP-even (CP-odd) Higgs bosons. As mentioned above 
we do not consider scenarios where the heavier H^ is produced in the cascades. The partial 
width for the neutralino decay (3.1) is given at tree- level by 

r(x° -> X°jH k ) = 16 ^3 rl/2 ( m I° , m|o , m 2 Hk ) (m|o + m|o - m 2 Hk + 2m^m^oj , (3.3) 

with a CP-even Higgs in the final state and 

\P "k\ 2 

for the decay into a CP-odd scalar. The Kallen function r(x, y, z) = (x — y — z) 2 — Ayz, 
and the coupling factors are 

e 



A 



(SkiN i3 - S k2 N i4 ) (c w N j2 - s w Nji) 



. ^5 (SttJV^ + S k2 N j3 ) + S k3 N H N j3 \ + -^=S k3 N i5 N j5 + 



(3.5) 
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and 



"Pijk 2cwsw 



(PkiN i3 - P k2 N i4 ) (c w N j2 - swNji, 



+ 

(3.6) 

— [iV i5 (P fcl iV i4 + Pk 2 N j3 ) + Pfes^iVys] - -^P k 3N i5 N j5 

where the mixing matrices Sij, Pij, and Nij are defined in section 2.1. Equations (3.3)- 
(3.6) assume a real neutralino mixing matrix and signed neutralino masses. Competing 
neutralino decay modes are into vector bosons, Xi — > x°j% an d Xi ~> xfW T . For brevity 
we refrain from giving expressions for these (and the corresponding chargino decay modes) 
here; they can be found in [57]. A detailed analysis of the mode is performed in [58]. 
Since the squarks and sleptons are assumed to be heavy, there are no open two-body decay 
modes of the neutralinos into the sfermion sector. Also slepton-mediated three-body decays 
- which can dominate over the two-body decays in certain scenarios — are numerically 
irrelevant for the same reason. 

The branching fractions for the relevant decay channels have been computed at leading 
order using FeynArts/FormCalc [59, 60] and a purpose-built Fortran code. 3 Results for the 
neutralino branching ratios in the modified P4 scenario are shown in figure 3. The decay 
modes of \\ and x% (upper row of figure 3) — which are both Higgsino-like — show similar 
patterns for large values of A K . The dominant mode is always Xi ~~ Xi% with a branching 
ratio of about 50%, but the Higgs channels are also significant with BR(x? — > Xi-^i) ^ 0-3 
and BR(Xi — > Xi-^2) ~ 0.15. An important point to note here is that the branching ratios 
of X2 an d X3 are quite insensitive to changes in Mh 1 (At)- For the heavier neutralinos, 
X4 and x® (lower row of figure 3), which also carry a larger gaugino fraction, the decay 
pattern is more complicated. Of largest interest for Higgs production is the sizable rate for 
Xa ~~ > X3H1 (once A K is sufficiently large to make this decay mode kinematically possible), 
and the fact that direct decays of X5 to the LSP are suppressed. This will lead to neutralino 
decay chains with intermediate (Higgsino) steps. Everything taken together, we can expect 
a large number of light Higgs bosons to be produced in neutralino cascade decays. 

The light chargino Xi decays exclusively into the LSP and a W boson, while the 
corresponding decay channels for the heavier chargino are shown in figure 4. Even if 
the dominant mode is X2 ~^ xfW^, independently of A K , there are several channels with 
a branching fraction of order 20% of interest for Higgs production. These include the mode 
X2 ~^ X1H2 a ud the decays into intermediate-mass Higgsinos, xf - > X2 s^ 1 * 1 - 

The heavier neutralinos and the heavy chargino, from which a Higgs could emerge as 
decay product, can either be produced at the LHC directly or in the decay of a heavier 
SUSY particle. The cross section for direct production of neutralino pairs is small, only 
O(fb) at y/s = 14 TeV, and the reach in these channels will be rather limited even for 
high luminosity. The large cross sections for production of strongly interacting sparticles 
(squarks and gluinos), on the other hand, are potentially more promising as a source of the 
heavier neutralino states and the heavier chargino. Exploiting cascade decays of this kind 



3 A FeynArts model file for the NMSSM has been obtained using FeynRules [61] and SARAH [62]. Details 
on this implementation will be presented elsewhere. 
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Figure 3. Branching ratios in the modified P4 scenario for x°i ~^ X^jZ (solid), x°l ~ * XjH\ 
(dashed), Xi — > (dotted), and Xi ~^ X°jA\ (dot-dashed). The color coding indicates the final 

state neutralino j = 1 (black), j = 2 (blue), j = 3 (magenta), j = 4 (green), or the chargino mode 
X° -> X^ T (red). 
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Figure 4. Branching ratios of xf m the modified P4 scenario. Neutralino final states xf ~~ ^ XiW^ 
(solid) are shown with the same color coding as for figure 3. Modes with a final state Xi are shown 
in red: xf ~> Xi z (solid), xf -> XiHi (dashed), xf ~> Xi H 2 (dotted), xf -> xf^i (dot-dashed). 



furthermore has the advantage that additional high-p^ jets are produced, which facilitates 
triggering and event selection. 

We use Prospino to calculate the NLO cross sections for production of pp — > gg, 
pp — > qq, pp — > qq, and pp — > gq according to [63], with CTEQ6 [64] parton distributions 
and a common choice of renormalization and factorization scales as the average mass of the 
final state (SUSY) particles. Numerical results are given in table 2 for the two centre-of- 
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Masses (GeV) 

Mg Mg 


o"lo (pb) 

gg qq gq qq s 


0NLO (pb) 

gg qq M qq £ 


= 7 TeV 


750 750 
1000 750 
1000 1000 


0.03 0.23 0.25 0.05 0.56 
0.002 0.19 0.06 0.05 0.31 
0.001 0.03 0.02 0.004 0.06 


0.07 0.27 0.39 0.08 0.82 
0.006 0.21 0.10 0.07 0.39 
0.005 0.04 0.04 0.006 0.08 


yjl = 14 TeV 


750 750 
1000 750 
1000 1000 
1500 1500 


1.18 1.67 5.20 1.06 9.11 
0.15 1.41 1.86 0.96 4.38 
0.14 0.42 0.87 0.18 1.61 
0.01 0.04 0.05 0.01 0.10 


2.21 2.06 6.78 1.53 12.6 
0.32 1.59 2.44 1.34 5.69 
0.31 0.51 1.19 0.26 2.27 
0.01 0.05 0.07 0.02 0.15 



Table 2. Total production cross sections for pp — > gg, pp — > qq, pp — > gq, and pp — > qq at LO 
and NLO SUSY-QCD from Prospino. The squark cross sections are summed over the four "light" 
squark flavours. No kinematic cuts have been applied here. 



mass energies 7 TeV and 14 TeV. The cross sections for pp — > tt and pp — > bb [65] are also 
calculated and included in the analysis, but since they turn out to be significantly smaller 
than a(pp — > qq) they are not shown in the table. In order to give some indication of the 
expected change in the number of events for different scenarios, the results are presented 
for several values of the squark masses, Mg, and the gluino mass, Mg. The mass ranges are 
selected to respect the published limits from ATLAS [66-70] and CMS [71-74] based on 
the 2010 data. Taking into account also the most recent results [75, 76], the Mg = 750 GeV 
case appears to be under some pressure. We present the results of our analysis below for 
the two cases Msusy = 750 GeV and Msusy = 1 TeV (the leading order squark masses 
are obtained from Msusy through eq. (2.12), to which higher order corrections are then 
added). 

The nearly mass-degenerate squarks decay preferentially into the SUSY-EW sector. 
Direct decays into Higgs bosons (or Higgsinos) are negligible for squarks of the first two 
generations due to the small Yukawa couplings. In contrast to the MSSM, the neutralinos 
also have a singlino component to which no squark couples. The left-handed squarks decay 
mainly into the wino, <jl — > W°q, q^ — > W q' , while the right-handed squarks decay mostly 
to the bino, q~R Bq. Numerically, this leads to squark decay modes listed in table 3 for 
the case with a soft scalar mass of Msusy = 750 GeV. The squark decay pattern for 
Msusy = 1 TeV is qualitatively similar. 4 Since the gaugino components are largest in the 
two heaviest neutralinos, the neutralinos produced in the squark decays tend to give rise 
to cascade decays with several steps. 

Finally, we note that the gluinos decay 'democratically' through g — >■ qq into all 
flavours, with rates governed only by the available phase space. 

The main numerical difference is an increase of BR(u L -»■ q'xf ) to 60% at the expense of a reduced 
BR(M£ -»• q'xf) = 3.9 x 1CT 2 . 
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Table 3. Branching ratios for the first and second generation squarks into neutralinos and charginos 
in the modified P4 scenario with .Msusy = 750 GeV. Results for channels with a branching ratio 
below 10 -4 are not shown. 

4 LHC analysis 

In order to assess whether the process discussed in the previous section can be useful as a 
Higgs search channel at the LHC we perform a Monte Carlo simulation. Here we use as 
benchmark the modified P4 scenario with the two different settings for the soft scalar mass: 
Msusy = 750 GeV and Msusy = 1 TeV. The DR value of the gluino mass parameter is 
set to M 3 = 1 TeV. We select A K such that M Hl ~ 40 GeV, which also affects M H2 , M Al 
and the branching ratios in the two cases as discussed in section 3. We have chosen this 
value of Mh 1 as an illustrative example of our scenario with 20 GeV < Mm < Mz and 
in order to make contact with the analyses of the "CPX hole" in the MSSM with complex 
parameters. Our results however depend only very mildly on the specific choice for Mn x - 
The simulation results are presented below both for LHC running at centre-of-mass energies 
of 7 TeV and 14 TeV. 

The squark and gluino-induced cascades in general give rise to a final state with high 
multiplicities and several hard jets, as well as large missing transverse momentum due to 
the presence of the LSP at the end of each decay chain. The minimal signal cascades 
(defined to be those with at least one Higgs boson present) generated by the production of 
a single squark or gluino correspond to 

q^qX°i -> qXiH k -> qx\bb, nj > 1, n b > 2, (14 a) 

g -)• gq -)• gqx® ->• 39Xi#fc -> gqX\bb, nj > 2, n b > 2. (14 b) 

Equations (14 a) and (14 b) show the minimum number of light and heavy flavour (b-) jets 
expected in the signal. Each event contains production of a pair of sparticles and their 
associated jets, meaning that the full signature for production of at least one H\ in the 
hadronic final state will be rij > 2, n b > 2. Since direct decays of the heavier (mainly 
gaugino) neutralinos into the singlino LSP are practically absent (cf. figure 3), most signal 
cascades will contain an intermediate Higgsino step which will add further particles in the 
final state. The typical jet multiplicity will also be higher due to additional QCD activity, 
in particular for gluon-initiated processes. 
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4.1 Event Generation 

For the event generation, we use MadGraph/MadEvent 4 . 4 . 44 [77] to calculate the leading 
order matrix elements for pp — > gg, qq, gq, qq, ti, and 66. The different event categories are 
weighted by the corresponding NLO cross sections to produce an inclusive SUSY sample. 
The resonance decay chains are then generated with PYTHIA 6.4 [78] using the NMSSM 
decay rates calculated above as input through the SUSY Les Houches accord [36]. The 
PYTHIA generator is also used to produce additional QCD radiation through initial- and final 
state parton showers, for parton fragmentation, and to generate multiple interactions for the 
underlying event. This produces fully dressed hadronic events which are passed through the 
fast simulation of the ATLAS detector performance implemented in the Delphes package 
[79]. 5 Hadronic jets are clustered using the anti-kx algorithm [80] with a jet radius measure 
of R = 0.4. 

Since for the lightest Higgs boson the decay to bb is favored, the probability i]b to 
correctly identify jets originating from bottom partons (6-tagging efficiency) becomes a 
crucial quantity for the analysis. Based on [81] we parametrize this efficiency as a constant 
% = 0.6 with respect to both the detector geometry and the jet energy scale. Only jets in 
the central tracking region |?7| < 2.5 can be tagged. The rate for misidentification as a 6-jet 
is assumed to be r] c = 0.1 for charm jets, and r\ q = 0.01 for jets produced by light quarks 
and gluons. The actual tagging algorithm implemented in the Delphes simulation is not 
based on a particular experimental method to identify 6-jets. The algorithm determines if 
a jet is close enough in AR to a "true" b parton. When this is the case, the efficiencies 
given above are applied to determine if the tagging is successful or not. 

4.2 Backgrounds 

Based on the event signature, SM production of ti with at least one hadronically decaying 
W boson (or additional jet activity) constitutes an irreducible background to the Higgs 
signal. We can a priori expect this to be the most important SM background since the 
scale for the SUSY-QCD processes is high (> 1 TeV). In principle there are other sources 
of background from production of W + jets (bb), Z+jets (bb), direct production of 66+ jets, 
or from QCD multijets. The cross sections for these processes are large compared to 
the signal cross section, with QCD multijets the largest and thereby potentially the most 
serious. However, for QCD jet production to constitute a background to the Higgs signal 
simultaneously a double misidentification of heavy flavour jets and a large mismeasurement 
of the missing transverse energy is required. It is furthermore difficult to simulate this 
background reliably, since extreme kinematical fluctuations — or experimental effects - 
would be necessary to produce the signal-like events. A detailed study of the experimental 
effects would require a full detector simulation, which is beyond the scope of the present 
paper. However, the dominance of the ti background over other SM processes, such as 
W + jets or Z + jets, for our final state has also been demonstrated experimentally by the 
results from SUSY searches with 6-jets and missing Et [68]. We therefore proceed under 

5 Running the same Delphes analysis with the "CMS" detector setup and similar parameters for jets and 
heavy flavour tagging, no significant differences are observed in the output. 
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Figure 5. Transverse momentum distribution of the hardest (non-6) jet at 7 TeV (left) and 14 TeV 
(right) for the inclusive SUSY sample with Msusy = 750 GeV(dashed) and SM ti background 
(solid). The histograms are normalized to unity. 

the assumption that the cuts devised to suppress the irreducible ti background will also be 
efficient for suppressing the other SM backgrounds as well. 

For the normalization of the ti background we use the NLO cross section a(pp —> 
ti) = 902 pb (y/a = 14 TeV) and cr(pp -»• ti) = 162 pb (y/i = 7 TeV), computed with 
the HATH0R package [82] for m t = 173.3 GeV and MSTW2008 PDFs [83]. In this way a 
consistent NLO normalization is used for both the signal and background events. The ti 
background is generated in the same Monte Carlo framework as already described for the 
signal. 

In addition to the SM backgrounds, the process we are interested in receives an im- 
portant background from the SUSY cascade itself. Any final state containing two 6-jets 
which do not result from an intermediate Higgs boson contributes to this background. 
Attempting to suppress the SUSY background events would require additional cuts that 
depend on the kinematics of the decay chains. This is something which may indeed be 
possible to devise once information on the supersymmetric spectrum has become available, 
but since we do not want to make any particular assumptions on the pattern of the SUSY 
spectrum, no selection will be applied aiming to reduce the SUSY background. Instead we 
will consider the inclusive bb mass spectrum directly after applying the cuts designed to 
reduce the SM background to determine if a Higgs signal can be extracted. 

4.3 Event Selection 

As a first step, we perform a preselection of the expected event topology, demanding rij > 2, 
rib > 2. All reconstructed jets are required to have a minimum px > 25 GeV. 

Figure 5 shows the pt distribution for the hardest jet in each event, comparing the 
inclusive SUSY events (with Msusy = 750 GeV) to the ti background. We show the results 
for the two cases y/s = 7 TeV (left) and y/s = 14 TeV (right). In order to illustrate the 
effect of applying cuts to this variable, each histogram is normalized to unity. From figure 5 
it is clear that the leading jet from the SUSY events has a much harder scale compared to 
the ti events. This can be understood as a result of the large boost obtained by the light 
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Figure 6. Transverse momentum distribution of the second hardest jet at 7 TeV (left) and 14 TeV 
(right) for the inclusive SUSY sample with Msusy = 750 GeV (dashed) and SM ti background 
(solid). The histograms are normalized to unity. 
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Figure 7. Missing transverse energy at 7 TeV (left) and 14 TeV (right) for the inclusive SUSY 
sample (dashed) and SM ti background (solid). The histograms are normalized to unity. 



quark jets originating from squark decays. It can also be seen that there is only a minor 
scaling difference in the jet pt distribution between the 7 TeV and 14 TeV cases. The 
same is true for the second hardest (light) jet, for which the corresponding pj- distribution 
is shown in figure 6. Similar differences between signal and background can be observed 
also for the third and fourth jet when they are present. 

With each cascade ending in the stable LSP, a large missing transverse energy 
is expected for the signal events. This distribution is displayed in figure 7, and shows 
indeed that the SUSY distribution peaks at high $ T values (> 200 GeV). This is therefore 
an important discriminating variable to suppress the background from ti events, where 
the missing transverse energy is due to neutrinos from leptonic W decays. As already 
mentioned, a hard cut on is a l so necessary to suppress the background from ordinary 
QCD multijet events and direct production of bb. A further advantage of the large $j> is 
that it can be used for triggering. 

The final kinematical distribution we are going to consider is displayed in figure 8. It 
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Figure 8. Distribution in AR(bb) at 7 TcV (left) and 14 TeV (right) for the inclusive SUSY 
sample (dashed) and SM ti background (solid). The histograms are normalized to unity. In events 
containing more than two 6-jets, all possible combinations have been included. 



Msusy = 750 GeV 


Inch SUSY 


Signal 


Background 


Generated events 


18 465 


8 261 


10 195 


Tij > 2, > 2 


4 048 


2 786 


1 262 


> 250 GeV, p£ t2 > 100 GeV 


2 436 


1 738 


698 


$ T > 150 GeV 


1 735 


1 211 


524 


min [AR(bb)} < 1.5 


1 014 


774 


240 


Total efficiency 


5.5 x 10" 2 


9.4 x 10" 2 


2.4 x 10" 2 



Msusy = 1 TeV 


Inch SUSY 


Signal 


Background 


Generated events 


20 671 


10 923 


9 748 


rij > 2, nj, > 2 


5 313 


4 344 


969 


> 250 GeV, p£ t2 > 100 GeV 


4 642 


3 828 


814 


$ T > 150 GeV 


3 705 


3 036 


669 


min [AR{bb)} < 1.5 


2 544 


2 170 


374 


Total efficiency 


0.12 


0.20 


3.8 x 10~ 2 



Table 4. Number of events remaining after each step of the event selection at yfs = 7 TeV. The 
SUSY events are classified as signal or background based on the presence of (at least one) Higgs 
boson in the decay chain. The total number of generated events in the inclusive sample is arbitrary. 



shows the separation in AR = (Ar/) 2 + (Acp) 2 between pairs of 6-jets. For events with 
rib > 2 all possible combinations have been included. The signal distribution is seen to peak 
near the minimum separation of AR = 0.4 set by the jet measure, while the ti background 
prefers the 6-jets to be more back-to-back and peaks at AR ~ ir. 

The precise cuts applied — and their effect on the event selection — are shown for 
the SUSY events in table 4 (for the 7 TeV case) and table 5 (14 TeV). Table 6 gives the 
corresponding information for the SM ti background. Note that the number of generated 
events in these tables does not correspond to any particular luminosity, but is rather 
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Msusy = 750 GeV 



Incl. SUSY Signal Background 



Generated events 


23 771 


10 874 


12 897 


rij > 2, nj, > 2 


5 009 


3 610 


1 399 


> 250 GeV, p*T 2 > 100 GeV 


3 287 


2 422 


865 


$ T > 200 GeV 


1 935 


1 400 


535 


min [AR(bb)] < 1.2 


991 


775 


216 


Total efficiency 


4.2 x 10~ 2 


7.1 x 10" 2 


1.7 x 10" 2 



Msusy = 1 TeV Incl. SUSY Signal Background 



Generated events 


20 232 


10 557 


9 675 


Tij > 2, > 2 


5 428 


4 338 


1 090 


> 250 GeV, p^ t2 > 100 GeV 


4 852 


3 924 


928 


$ T > 200 GeV 


3 392 


2 719 


673 


min [AR(bb)] < 1.2 


1 983 


1 673 


310 


Total efficiency 


9.8 x 10~ 2 


0.16 


3.2 x 10~ 2 



Table 5. Events remaining after each step of the event selection at ^/s = 14 TeV. The event 
categories are similar to those in table 4. The total number of generated events in the inclusive 
sample is arbitrary. 

selected to give adequate statistics for the event selection. The inclusive SUSY sample is 
split into signal and background, where the signal consists of the events containing at least 
one Higgs boson (as determined from Monte Carlo truth information). In the last row we 
give the accumulated total efficiencies of all the cuts. Looking first at table 4, we see that 
an efficiency of 5.5 x 10~ 2 is obtained for the case with Mstjsy = 750 GeV. This efficiency 
is more than doubled (0.12) for the case with Mstjsy — 1 TeV, since the heavier squarks 
give harder jets as decay products which leads to more events passing the jet pr cuts. The 
larger boost given to the LSP at the end of the decay chain also leads to an increased pj>. 
The same qualitative features are visible at 14 TeV, as can be read off table 5. Due to the 
favorable signal statistics at 14 TeV, 6 we can afford slightly harder cuts on $j> and AR(bb) 
in this case, something which is also needed to maintain a good background suppression. 
One should therefore not be discouraged by the somewhat lower efficiencies recorded in 
this case (4.2 x 10~ 2 for M SU sy = 750 GeV vs. 9.8 x 10" 2 for M SU sy = 1 TeV). The signal 
efficiencies can be compared to those for the ti background, given in table 6, which are 
at the 10 -5 level for both energies. It is clear from this table that the hard cuts on the 
jet pt and the distribution are the most important handles available to suppress the 
background. 

As discussed in the previous section, we do not apply any specific cuts to suppress 
the background from SUSY events that do not involve a Higgs boson. The numbers given 
in tables 4 and 5 show that nevertheless our event selection gives rise to an improvement 
also in the ratio of signal events over SUSY-background events. The largest difference in 

"Going from 7 TeV to 14 TeV, the signal cross section for Msusy = 750 GeV increases by a factor 14.5, 
while the ti cross section is only increased by a factor 5. 
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SM tt background 



7 TeV 14 TeV 



Generated events 
rij > 2, rib > 2 

> 250 GeV, ^ t2 > 100 GeV 
$ T > 150 GeV (7 TeV), > 200 GeV (14 TeV) 
mm[AR(bb)] < 1.5 (7 TeV), < 1.2 (14 TeV) 



900 000 



259 110 



1 120 



102 



12 61 



2 000 000 
576 232 



5 189 
405 



Total efficiency 



1.3 x 10" 5 3.0 x 10" 5 



Table 6. SM tt background events remaining after each step of the event selection at y/s = 7 TeV 
and yfs = 14 TeV. The total number of generated events is arbitrary. 

selection efficiency between the SUSY signal and background arises from the typical number 
of b quarks produced in the two cases, which is larger for the events where Higgs bosons are 
produced, leading to a stronger reduction of the SUSY background by the jet multiplicity 
cut. The cut on AR also contributes to the difference. This cut has the pleasant "side 
effect" to enrich the SUSY sample in Higgs events since the jets resulting from H\ — > bb 
decays are more likely to show up for small AR than those from two unpaired 6-jets. 

4.4 Results 

Figure 9 shows the resulting bb mass spectra after final event selection for an integrated lu- 
minosity of 5 fb" 1 at 7 TeV. For events with n\ ) > 1 only the fr-jet combination minimizing 
AR(bb) has been included. This reduces effects of combinatorics and increases the sensi- 
tivity for discovering resonances in the low mass region. For the scenario with relatively 
light squarks (Mstjsy = 750 GeV) shown in the left plot, we observe two peaks close to the 
masses of the Z boson and Hi, respectively. There is also a continuous distribution with 
a tail towards much higher values for M b g. This results from false pairings, fake 6-jets, or 
from b-jet pairs of non-resonant origin such as t or b decays. The same qualitative features 
are visible in the signal for Mstjsy = 1 TeV (right plot), but the statistics is rather poor 
due to the low signal cross section. In figure 10 we show the M b i distribution at 14 TeV, 
again for an integrated luminosity of 5 fb _1 . Here the signal statistics is much higher, so 
that a clear distinction of the H% resonance from the background should be possible both 
for Mstjsy = 750 GeV (left) and M SU sy = 1 TeV (right). 

The same distributions are shown in figures 11 (for the LHC at 7 TeV) and 12 (for 
14 TeV), but here with stacked histograms to more closely resemble "real" data. Here we 
have furthermore split up the inclusive SUSY sample into signal events (displayed in red), 
characterised by the presence of (at least) one Higgs boson in the decay chain, and the 
remaining SUSY background events (black). The latter constitutes an additional source of 
background besides the SM ti background (light gray). In figure 11 we see that the most 
striking feature is the H\ peak. Although the ti background peaks at roughly the same 
position as the signal, the statistics of signal events should be sufficient for establishing a 
signal over the background. In the 14 TeV case, figure 12 illustrates the features observed 
already in figure 10. The H\ peak stands out clearly above the background distribution, 
both for Msusy = 750 GeV and M SU sy = 1 TeV. 
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Figure 9. Invariant mass of 6-jet pairs after final event selection for the inclusive SUSY sample 
(red, dashed) and SM tt background (solid) at 7 TeV in the modified P4 scenario with Msusy = 
750 GeV (left) and Mgusy = 1 TeV (right). The histograms have been normalized to an integrated 
luminosity of 5 fb _1 . 
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Figure 10. Invariant mass of 6-jet pairs after final event selection for the inclusive SUSY sample 
(red, dashed) and SM tt background (solid) at 14 TeV in the modified P4 scenario with -Msusy = 
750 GeV (left) and Msusy = 1 TeV (right). The histograms have been normalized to an integrated 
luminosity of 5 fb _1 . 



With the assumed statistics, no peaks are observed in the bb mass spectrum for the 
heavier Higgs bosons H2 and A\. This is mainly due to the smallness of the branching 
ratios into the bb mode because of the open Higgs decay channels. Part of the difficulty in 
observing the heavier resonances is also a result of selecting the combination minimizing 
AR(bb) in configurations with multiple fo-jets, which favors selection of the light H±. 

In order to obtain an estimate of the significance of the H± mass peak we have per- 
formed a Gaussian fit to the maximum of the distributions in figures 11 and 12. Table 7 
lists the results extracted from the fit for the mean value Mjj and the 1 a width AMh of the 
Gaussian peak. We find that the fitted central values reproduce well the correct Hi mass 
for all cases (recall that the input mass used in our numerical simulation is Mjj 1 = 40 GeV). 
The statistical uncertainty on the mean value Mh from the fit is about ±3 GeV for the 
LHC energy of 7 TeV and about ±1 GeV for 14 TeV. This reflects the lower signal statis- 
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Figure 11. Invariant mass of 6-jet pairs for SUSY signal (red), SUSY background (black) and 
SM tt background (light gray) in the modified P4 scenario with Msusy = 750 GeV (left) and 
-Msusy = 1 TeV (right) at 7 TeV for an integrated luminosity of 5 fb _1 . 
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Figure 12. Invariant mass of 6-jet pairs for SUSY signal (red), SUSY background (black) and 
SM tt background (light gray) in the modified P4 scenario with Msusy = 750 GeV (left) and 
-Msusy = 1 TeV (right) at 14 TeV for an integrated luminosity of 5 fb -1 . 



tics available in the low energy running, and the more coarse binning in M bb required to 
observe the peak. 

The number of signal and background events in the peak region is obtained by inte- 
grating Mbb over the interval [Mjj — AMjj, Mjj + AMjg-], corresponding to drier of the 
Gaussian distribution. As explained above, the combined background includes both the 
events from SM tt and the part of the inclusive SUSY sample containing no Higgs bosons 
in the cascades. The event numbers are combined into the ratios of signal/background 
(S/B) and S/VB given in table 7. We use S/\fB as a simple illustration for the expected 
significance and in particular for comparing between the four example cases we consider 
here and with other theoretical studies using the same criterion. Clearly, claiming an 
actual discovery would require a more sophisticated statistical treatment. We regard it 
nevertheless as encouraging that a significance of S/y/B > 5 is achieved for three of the 
four cases considered in table 7. The only exception is the case Msusy = 1 TeV for the 
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M SUSY (GeV) M H (GeV) AM H (GeV) S/B S/VB 



7 TeV 
7 TeV 
14 TeV 
14 TeV 



750 41.6 12.2 3.4 8.2 

1000 37.7 17.6 1.4 3.9 

750 39.5 8.0 4.5 29.7 

1000 39.4 9.7 4.3 29.3 



Table 7. Mean value Mjj and width AMjj extracted from a Gaussian fit to the bb mass peak. The 
corresponding number of signal (S) and background (B) events is recorded in the ±AMh interval 
around the fitted resonance mass. All results are presented for an integrated luminosity of 5 fb _1 . 

LHC at 7 TeV, where S/\Z~B = 3.9 with 5 fb _1 . A somewhat higher luminosity (or a 
combination of ATLAS and CMS data) would be needed in this case in order to reach a 
significance S/y/B > 5. The very high significances of about 30 obtained for the 14 TeV 
case illustrate the qualitative features already observed in the discussion of figures 10 and 
12: there should be no problems in establishing a signal in this mass region. 

5 Summary and Conclusions 

The NMSSM is both theoretically appealing as an extension of the SM and interesting 
phenomenologically, as its spectrum may contain Higgs bosons with mass much below the 
limits in the SM or the MSSM. We have investigated an NMSSM scenario with a light 
CP-even Higgs in the mass range 20 GeV < Mn x < Mz- Scenarios like this may be missed 
with the standard Higgs search channels at the LHC, in particular due to a potentially large 
branching ratio of the heavier H2 state, that has SM-like couplings to gauge bosons, into 
a pair of light Higgses. We have pointed out that there are good prospects for discovering 
such a light Higgs boson in SUSY cascade decays at the LHC. 

We have performed a Monte Carlo simulation of the signal and the dominant back- 
ground to the level of fast detector simulation, taking into account also background from 
other SUSY events that do not involve cascade decays containing a Higgs boson. For 
our numerical analysis we adapted the "P4" benchmark point proposed for the NMSSM, 
choosing Mg x = 40 GeV as example value for the mass of the light Higgs. Production of 
squarks and gluinos via the strong interaction at the LHC may give rise to cascade decays 
involving heavy neutralinos and charginos decaying into lighter ones and a light Higgs. We 
have investigated the impact of various kinematical variables on discriminating between 
the inclusive SUSY signal (including events both with and without a Higgs boson in the 
cascade) and the SM background from ti production. A set of simple cuts has been devised 
that turned out to be efficient for establishing the inclusive SUSY signal. We did not as- 
sume any specific knowledge about the background from SUSY events without a Higgs in 
the cascades. Accordingly, besides favoring events containing the light Hi by selecting the 
combination minimizing AR(bb) in configurations with multiple 6-jets, we have not applied 
any particular cuts for suppressing the SUSY background. 

Our results show that reconstruction of the decay of the light Higgs into bb may be 
feasible. Such an observation would be a direct experimental sign of the bottom Yukawa 
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coupling, which is difficult to access in standard search channels. We have investigated two 
values of the soft SUSY-breaking parameter in the squark sector, Msusy = 750 GeV and 
MgusY = 1 TeV, while we set the gluino mass parameter to 1 TeV. A modest integrated 
luminosity of 5 fb _1 has been considered for LHC running both at 7 TeV and 14 TeV. We 
find a statistical significance for the H\ mass peak of S/V~B ~ 4 for Msusy = 1 TeV at 
\fs = 7 TeV. This significance increases to S/VB ps 8 for A/susy = 750 GeV at 7 TeV 
and reaches a level of almost 30 for both values of Msusy at 14 TeV. While the example 
values that we have chosen for Msusy an d the gluino mass are close to the current search 
limits from the LHC, the large statistical significance that we have found for the 14 TeV 
case indicates that there is certainly scope to extend our analysis to scenarios with heavier 
squarks and gluinos or to scenarios with reduced branching ratios of the neutralinos into 
Higgs bosons. Since the high-energy run of the LHC is not imminent, we leave a more 
detailed analysis of this reach for future work. 

The results presented here have been obtained in a specific benchmark scenario, but it is 
easy to see that they are more generally applicable. First of all, the value Mh 1 = 40 GeV 
used in our numerical analysis was chosen just for illustration. Our results are rather 
insensitive to the precise value of Mh\ ■ Since the production relies on the decay of heavier 
SUSY states, with branching ratios largely independent of M/^, the Higgs production 
rates remain similar for the whole mass range Mjj-y < Mz- The event selection and signal 
identification through H\ — > bb proceeds along similar lines as we have discussed. 

Concerning the settings of the other SUSY parameters, our results will be similar for 
other scenarios fulfilling a few simple criteria: Obviously, the neutralinos and charginos 
have to be sufficiently lighter than the squarks and gluinos in order to be produced at all 
in the cascade decays of the latter. The squark decays also provide the hard jets utilized 
in the event selection. With the present limits from the LHC searches on the masses of the 
gluino and the squarks of the first two generations this criterion is almost automatically 
fulfilled for any model of interest. Furthermore, the neutralino and chargino mass hierarchy 
and mixing character must be such that the squark decays go through heavier neutralinos 
or charginos, and the decays of the latter into a light Higgs and a lighter neutralino or 
chargino are open. Such a scenario is disfavored if the LSP is gaugino-like. In order to 
generate a sufficient number of Higgs bosons in the cascade decays, it is also advantageous 
for (at least one of) the gauginos to be heavier than the Higgsinos, so that an intermediate 
Higgsino decay step can be present. In the NMSSM such a situation can be realized quite 
easily if the LSP is singlino-like. 

While the results presented in this paper are based on a rather simple-minded analysis, 
involving for instance just a fast detector simulation, we nevertheless regard them as very 
encouraging, motivating a further exploration of the potential for detecting a light non-SM 
type Higgs in SUSY cascade decays. In fact, there exists the exciting possibility that the 
discovery of a SUSY signal could go hand in hand with the discovery of one or more Higgs 
bosons. 
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